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Abstract. Brownian motion is a well-known model for normal diffusion, but not all physical phenomena behave accord- 
ing to a Brownian motion. Many phenomena exhibit irregular diffusive behavior, called anomalous diffusion. Examples of 
anomalous diffusion have been observed in physics, hydrology, biology, and finance, among many other fields. Continuous- 
time random walks (CTRWs), introduced by Montroll and Weiss, serve as models for anomalous diffusion. CTRWs gen- 
eralize the usual random walk model by allowing random waiting times between successive random jumps. Under certain 
conditions on the jumps and waiting times, scaled CTRWs can be shown to converge in distribution to a limit process M{t) in 
the cadlag space ID)[0, <») with the Skorohod J[ orM| topology. An interesting question is whether stochastic integrals driven 
by the scaled CTRWs X"{t) converge in distribution to a stochastic integral driven by the CTRW limit process M{t). We 
prove weak convergence of the stochastic integrals driven by CTRWs for certain classes of CTRWs, when the CTRW limit 
process is an a-stable Levy motion and when the CTRW limit process is a time-changed Brownian motion. 

Key words: Continuous-time random walks. Stochastic integrals. Weak convergence. Anomalous diffusion, a-stable Levy 
motion. Time-changed Brownian motion. 



1. Introduction 

Continuous-time random walks (CTRWs) were developed by Montroll and Weiss fTSl to study random walks on a 
lattice. CTRWs generalize the usual random walk model by introducing random waiting times between jumps. The 
addition of random waiting times allows CTRWs to capture irregular diffusive behavior, making CTRWs good models 
for anomalous diffusion. Unlike normal diffusion, anomalous diffusion is characterized by a non-linear relationship 
between the mean squared displacement of a particle and time. Anomalous diffusion has been observed in hydrology, 
where contaminants often travel more slowly in ground water due to sticking or trapping, in biology, where proteins 
diffuse more slowly across cell membranes, and in the behavior of many financial markets where heavy-tailed price 
jumps occur lfT0l[T6l . Additional examples of anomalous diffusion can be found in physics, hydrology, biology, and 
finance, among many other fields flOl [161 [iTl l20l l22l . Because of the existence of so many examples of anomalous 
diffusion, it is important to be able to model and to understand this type of behavior 

Physicists often use the CTRW model to derive a partial differential equation which describes the anomalous dif- 
fusion. The differential operators in the normal diffusion equation ^^^7^ = D^-j^^ may be replaced by fractional 
derivatives in time and space or even pseudo-differential operators in time and in space lfT5l[T6ll . The density of the 
"long-time scaling limit" of the CTRW solves the space-time-fractional diffusion equation derived from the CTRW 
model ifTSl flTl. Moreover, the long-time scaling limit of the CTRW is actually the weak limit of a sequence of normal- 
ized CTRWs. Meerschaert and Scheffler derive the scaling limit of the CTRW under certain conditions and identify 
the CTRW limit process as a time-changed Levy motion (see fTsl for details). Studying the CTRW limit process 
and the processes it drives is important, since the dynamics of the limit process correspond to its anomalous diffusion 
equation. 

Many phenomena are modeled by a stochastic differential equation driven by the CTRW limit process (for example, 
see lU, 14 1). For applications, it is also important to consider how the stochastic integrals driven by the CTRW limit 
process can be approximated. Let Bj^* denote the cadlag function space of M'^-valued functions, J\ the standard 
Skorohod topology on D |2|, and => convergence in distribution. We study the following question: 
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Q: lfX"{t) is a scaled CTRW and if {H",X") {H,M) in {B^2,Ji) where M{t) is the CTRW limit process, then 
under what conditions does 

(*) \ f H"is-)dX"{s)\ f H{s-)dM{s)\ in(DK,7i)asn^-? 

A paper of Germano et al. conjectures that for the usual compound Poisson process X{t), the stochastic integrals 
of the scaled process X"{t) driven by X"{t) converge weakly to the stochastic integral of Brownian motion driven by 
Brownian motion [31, i.e. {/^^"(i— )c/X"(i)},>o =^ {fQB{s)dB{s)}t>a- In this paper, we verify this conjecture as a 
special case of a more general result (see Example|2|. We show that © holds in the case that the limit process is an 
a-stable Levy motion (Theorem |2]l, and in the case in which the limit process is a time-changed Brownian motion 
(Theorem[3]). In order to prove these results we rely on Kurtz and Protter's theorem on weak convergence of stochastic 
integrals [ 13 1. The key condition is that the integrators {X") be uniformly tight. 

In Section m we recall the required facts on continuous-time random walks, scaling limits, and weak convergence 
of stochastic integrals. Our main results are contained in Section[3]as well as the example proving the conjecture of 
Germano et al. Section|4]focuses on some applications and examples of the main results. 



2. Background 

2.1. Continuous-time random walks. Under the random walk model, random displacements occur at regular deter- 
ministic intervals. The continuous-time random walk (CTRW) model allows random waiting times between successive 
random displacements, or jumps. Let (^,),gn be a sequence of independent and identically distributed (i.i.d.) M-valued 
random variables representing the particle jumps and (t, be a sequence of i.i.d. nonnegative random variables rep- 
resenting the waiting times between successive jumps. Throughout, the focus is on uncoupled CTRWs in which the 
sequences of jumps ((^,) and waiting times (t,) are independent. Set 5(0) = and r(0) = and let S{n) = + . . . + ^„ 
be the position of the particle after « jumps and T {n) = T\ + . . . + T„ be the time of the nth jump. For f > 0, define 

(1) A^(f) :=max{« >0: r(«) <f} 

to be the number of jumps by time f . The position of the particle at time t can then be written as the sum of the jumps 
up to time f : 

m 

(2) X{t)^S{N{t))^Y^^t. 

1=1 

The stochastic process {X(f)},>o is called a continuous-time random walk (CTRW). 
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Figure 1 . A continuous-time random walk. 



Since the waiting times (t,) are i.i.d., N{t) is a renewal counting process and the CTRW is a compound renewal 
process, or renewal-reward process (as seen in queueing theory). A CTRW may also be called a point process with 
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reward. The counting process N{t) and the the time of the nth jump T{n) have a particularly nice inverse relationship. 
For any n>Q and any f > 0, 

(3) {N{t)>n} = {Tin)<t}. 

Equation (O holds because the number of jumps by time t is at least n if and only if the nth jump occurs at or before 
time t. 

Every CTRW is a semi-Markov process and is Markovian if and only if the distribution of the waiting times is 
exponential 11211 1711. Germano et al. |21 1 prove that a CTRW is also a martingale under certain conditions. Since we 
require this result in Section[3] we include an alternate (and more direct) proof here. 

Lemma 1. Let X{t) = Y!!=Hi be a CTRW. Assume E{E,i) =0, E\E,i\ < oo, and E[N{t)] < oo far each t. ThenX{t)is 
a martingale. 

Proof. First observe that X(f ) € L' for all f > by Wald's Equality, 



!=1 



N(t) 

<EY^\^,\^EME[N{t)]« 

i=\ 



Now let = (y{X{s),Q <s< f}, the natural filtration of X. ThenX(f) is adapted to and for s < f, 

E[X{t)\.^,;\=E[X{t)-X{s)\.^,]+E[X{s)\ 

N{t) 

_i=N[s)+l 
- Nit) 1 

_i=N(s) + l 

'-E[^i]{E[N{t)] -E[N{s)]) +X{s) by Wald's Equality 
--X{s) since £[(^1] =0. 



Xis) 

X{s) since ^; are i.i.d. 



□ 



lfX{t) is a CTRW satisfying the conditions in Lemma[T]and//(f) a cadlag adapted process, thenX(f) is a martingale 
and we can understand jQH{s—)dXs in the sense of Protter 1 19|. But we can also make a simpler definition. Germano 
et al. [21 1 note that this integral can be defined as a sum, since the CTRW is a step function, i.e.. 



(4) 



H{s-)dXis) 



N{t) 

!=1 



where the integrand is evaluated at the jump times of X. The right-hand side can be considered as a coupled CTRW, 
since the new jumps H{Ti—)^i clearly depend on N{t). 

2.2. Scaling limits. Here we recall the known scaling limit theorems for CTRWs. By the scaling limit of a CTRW, 
we mean the limit process resulting from appropriate scaling in time and space according to a functional central limit 
theorem (FCLT). The limit behavior of the CTRW depends on the distribution of the jumps and the waiting times. If 
the waiting times have finite mean, the CTRW behaves like a random walk in the limit. So, by Donsker's Theorem, if 
the waiting times have finite mean and the jumps have finite variance then the scaled CTRW converges in distribution 
to a Brownian motion. If the waiting times have finite mean and the jumps are in the domain of attraction (DOA) 
of an a-stable random variable, with a e (0,2), then the appropriately scaled CTRW converges in distribution to an 
a-stable Levy motion ll23l Theorem 4.5.3]. 

However, if the waiting times have an infinite mean, the CTRW limit behavior is more complex. For this case, 
Meerschaert and Scheffler prove a FCLT which identifies the limit process as a composition of an a-stable Levy 
motion A{t) and the inverse of a j3-stable subordinator E{t), where a E (0,2] and j3 G (0, 1) (see Theorem 4.2 in 
[TSl for details). Their proof uses a continuous-mapping approach and the convergence holds in the Mi -topology 
and not in the stronger -topology. By examining the proof, we note that their convergence result only holds in the 
Ml -topology since the composition map is Mi but not 7i -continuous at{A,E). However, in the case that a = 2, we can 
prove that the the convergence holds in the J\ -topology. This is important since the theorems on weak convergence of 
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stochastic integrals are stated in the Ji -topology. We provide a proof (in W) for the case a — 2 here following the proof 
of Meerschaert and Scheffler The resulting limit process, a time-changed Brownian motion B{E,), is of particular 
interest since it is a good model for subdiffusion. 

Theorem 1. Let (t,),^^ be i.i.d. nonnegative random variables in the strict domain of attraction of a ^-stable, 
j3 G (0, 1). Let (^,),gn be i.i.d. R-valued random variables in the strict domain of attraction of a normal (i.e., an a- 

^' 

stable with a — 2). Let X{t) = ^ ^j. Then there exists a regularly varying function c of index j3 /2 with c{n) — c(b{n)), 

i=\ 

b regularly varying of index j3, and c{n) — ^JnL{n) for some slowly varying function L such that, as n ^ °°, 

^{B(£(0)}r>0 /«(BM,/i). 



Proof. We follow the continuous-mapping approach used in ||T5l . substituting the 7i -continuity of composition the- 
orem instead of the Mi. Let B"{t) = ^ } and E"{t) = . By the independence of the jumps (^,) and the 

c(n) bin) 

waiting times (t,). Theorem 4. 1 and Corollary 3.4 in ifTSl imply that as « ^ oo, 

{(B"(f),£"(f))},>o ^ {B(0,£(0}r>o in (Pr x Dr^,7i). 

Consider the continuous-mapping theorem fTi), Theorem 3.4.3] where the map is composition, o : (Dr x — > 
(Dr,7i). Let C be the space of continuous functions, D| the space of nondecreasing cadlag functions, and C|| the 
space of strictly increasing continuous functions. By the 7i -continuity of composition theorem ll23l Theorem 13.2.2], 
the set of discontinuity points of composition under Jy is Disc(o) = ((C x D|) U (B x Cpf))'^. Since {B{t),E{t)) G 
C X D|, P{{B,E) e Disc(o)) = 0. Since o is measurable by |2, p. 145], the continuous-mapping theorem implies that 
B" o E" ^ B o E in the Ji topology. Consequently, as n — > oo, 

'X{nt)' 



c{n) 



^ {B(£(f))}r>o in(BM,7i^ 

t>0 



□ 



2.3. Weak convergence of stochastic integrals. We now review the conditions under which a sequence of stochastic 
integrals converges in distribution. Let 0" — (^"),>o,-P") be a sequence of filtered probability spaces satis- 

fying the usual hypotheses: (1) contains all the P"-null sets of and (2) the filtration (^") is right-continuous, 
i.e., = f]e>o-^"+£- Assume (//") is a sequence of cadlag ^"-adapted processes, X" a sequence of cadlag -J^,"- 
semimartingales, and X" X. Kurtz and Protter call {X") good if for any sequence (//") in Dr defined on 0" such 
that {H",X") ^ {H,X) in (Djj2,/i), then there exists a filtration {^t) such that X is an (^()-semimartingale and 
/ H'^_ dX^ / dXs in (Dr,7i) |12, Definition 13]. 

Necessary and sufficient conditions for goodness are discussed in |T3lfT2ll . First, it is necessary to assume [H" ,X") 
converges jointly to {H,X) in Djj2 and not in Dr x Dr. If this joint convergence holds, then the key condition is that 
the integrators {X") be uniformly tight (UT) or, equivalently, have uniformly controlled variations (UCV) 111 21 . In fact, 
if {H",X") ^ [H,X) in {D^2,Ji), then {X") is good if and only if (X") satisfies UT (equivalently UCV) lEl Theorem 
32, 34]. We use the UT condition throughout this paper The definition of uniform tightness was first introduced by 
Jakubowski, Memin, and Pages |9|. A sequence of semimartingales iX") is said to be uniformly tight (UT), if for each 
t > 0, the set { jl,H"_dX" ,H" e S'Mi/"| < \ ,n> 1} is stochastically bounded (uniformly in n), where S" denotes the 
collection of simple predictable processes on 0". Since we desire the weak convergence of the stochastic integrals 
/ H"_ dX" f Hs- dXs, it makes sense to require tightness of the laws of / H"_ dX". However, it is not easy to verify 
the UT condition directly. We rely on the following proposition, which gives a more accessible condition for UT to 
hold. Let |AX"(f)| = |X"(f) -X"(f-)| be the size of a jump at time f. 

Proposition 1. Il9l ll2l If (X")„yi is a sequence of local martingales and if for each t <°°, 



(5) sup£" 

n 

then {X") satisfies UT. 



sup I AX" (i) I 
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3. Results 



Our first result concerns the weak convergence of a stochastic integral driven by a CTRW with deterministic waiting 
times to a stochastic integral driven by the CTRW limit process A (f), an a-stable Levy motion. Let S'a = 5a(l, 7,0) be 
the family of a-stable distributions with skewness parameter ye [—1,1]. Our proof of the following theorem requires 
that the jumps have a finite first moment, so we assume a E (1,2]. This assumption is sufficient for most applications 
involving superdiffusion. 

Theorem 2. Let be i.i.d. mean random variables. Assume belongs to the domain of attraction of an 

a-stable random variable Sa, OC G (1,2]. Define 

a(n) ain) 

where a{n) ~ n^^"L{n) with L a slowly varying function. If(H",X") {H,A) in (Dg2,/i), then there exists a filtration 
(■^t) such that A is an {■^t)-semimartingale, H is an {^t)-odapted cadldg process, and as n °° 

(7) \ f H"{s-)dX"{s)\ f H{s^)dA{s)\ /«(©R,7i). 



Our goal is to show Theorem|2]follows from Proposition[T]and Theorem 34 in lfT2l once the conditions of Proposi- 
tion[T]are verified. 

Proof. We first show that {X"{t))„>i is a sequence of martingales. Note that the jumps i*, have mean and finite first 
moment since a > 1. By the same techniques used in Lemma[T] we see thatX"(f) is an ^"-martingale (and therefore 
a local martingale) for each n, where is the natural filtration of X". Therefore, the first hypothesis of Proposition 
□holds. 

To verify the second hypothesis, we need to show that sup,, £" [sup^.<, |AX"(s)|] < oo for each t < oa. Fix f < oo. 
Observe that, since X"(s) is a cadlag step process, it has jumps of size \ -^ \ at times k/n, A: G N. Additionally, since 
X" (f ) has finitely many jumps by time f, we can replace the sup up to time f by a max up to \nt\ . So we need to find for 
each / a uniform bound in n for E" maxj<^.< | | . Let (e^.) be a sequence of i.i.d. Rademacher random variables, 

[nt\ 

independent of {^k)- Then |<^<.| = \£k^k\ ™d Ej^^i^ are i.i.d. in DOA{Sa)- Let S{nt) = ^ e^^^, the symmetrized sum. 

k=\ 

By Levy's Inequality, 



E" 





^k 






\ 


£k^k 




max 


a{n) 




max 


a{n) 


\<k<[nt\ 






l<*:<[n(J 





< I 2P"< 



k^i fl(n) 



> X = 2£'" 



S{nt) 



a{n) 



Note that for each n. 



2E" 



S{nt) 



a{n) 



[nr\ 

< 2 ^ £" 

(=1 



a{n) 



2[nt\E" 



a{n) 



<^ oo. 



^1 has a finite first moment since it is in the domain of attraction of Sa, CH E (1,2]. So, for n small, the expectation is 
bounded. It remains to show the expectation is bounded in the case of large n. Observe that 



2£" 



S{nt) 



a{n) 



Lin) 



S{nt) 



a{\nt\ ) 



Since e^^^, are in DOA{Sa) for a e (1,2], then by |[T1 Exercise 9, p. 91] the following first moments converge: 



(8) 



S{nt) 
<\nt\) 



■E\Sa\ as n — ^ oo. 
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By applying Equation (O together with a corollary by Feller in (A^.. p. 274] on slowly varying functions, there exists 
some «o G N such that for all n > no. 



2E" 



S{nt) 



a{n) 



<2cfi/«(£|5„| + l)<- 



for c > 1 , where Sa has a finite first moment since a G (1,2]. 
So the desired inequality holds: 



sup£"i sup|AX"(i)| I =max( max £"( sup |AX"(5)| |, sup £"( sup |AX"(5)|| 



< max 2[nt\E" 

l<n<nQ 



a{n) 



-2rti/«(£|5„| + l)<. 



Since {X") satisfies UT, we apply Theorem 34 in lE) to conclude that if {H",X") ^ {H,A) in (Bri then there 
exists an appropriate filtration such that as n — > oo, 

fH"{s~) dX'\s)\ =^ ( f H{s-) dA{s)\ in (Br,7i). 

□ 

By an easy application of this theorem, we see that the stochastic integral of a scaled simple random walk with 
respect to itself converges in distribution to the stochastic integral of Brownian motion with respect to Brownian 
motion. 



\nt\ 



Example 1. LetX"(f) = where t,i are i.i.d. mean random variables with variance 1 . By Donsker's Theorem, 

{^"(f)}r>o ^ {B{t)}t>o in (Bm,/i), so {X'\X") ^ {B,B) in (©1.2,71). Then by Theorem|2]as n ^ 00, 



(9) 



fx"{s-) dr\s)\ =^ ( f B(s) dB{s)\ in (Bffi,7i). 

^0 J r>o Uo J ,>o 



We now establish a result concerning the weak convergence of stochastic integrals driven by CTRWs to a stochastic 
integral driven by the CTRW limit process B{Et) when the CTRWs have infinite mean waiting times and jumps in the 
domain of normal attraction of a normal law. Specifically: 

Theorem 3. Let ((^,),gn be i.i.d. mean random variables with E\E,\ p = c^, so t,i belongs to the domain of normal 
attraction of a normal randomvariableZ^N{Q^\). Let (Tj)^^^^ be i.i.d. strictly P -stable random variables, fi E (0,1). 
Define 



(10) 



Nm p. 

x«(f)-E 

i=l 



/3/2- 



If{H",X") {H,B{E)) in (B^2,7i), then there exists a filtration {^t) such that B(E) is an {j^t)-semimartingale, H 
is a {.^t)-odapted cddldg process, and as « — > 0° 



(11) 



H"{s-)dX"{s) 



r>0 



H{s-) dB{E,) 



in (Bffi,7i) 



f>0 



As in the proof of Theorem |2] we prove this result by first verifying that the conditions of Proposition [T] hold and 
then applying Theorem 34 in {Vl\. However, since this case involves random waiting times with infinite mean, our 
proof also requires a uniform bound on the expectation of the scaled counting process N„t jn^ . We prove the necessary 
lemma here before presenting the proof of Theorem[3] 

Lemma 2. Let (t,),gn be i.i.d. strictly ^-stable random variables, j3 £ (0, 1). Then for each t >Q, 



(12) 



for all n> I, where M is a constant. 
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Proof. The key to bounding the expectation will be the asymptotics of the density of Ti at zero. To get at Ti , we write 
the expectation as the integral of the tail probabilities of the process A',,, / rfi . Using the inverse relationship given by 
Equation (O between A'r and T{n) as well as the fact that the T, are strictly j3-stable, we obtain 



P{N„, > xn^ ) dx 
p{T{\xnPA <nt] dx 



nt 



dx. 



Since mx'/^ < [xn^^] , the expectation is bounded above by P < -^rjp^ dx. By change of variables, this 
integral can be re-written in a more useful form. Explicitly: 



dx ~ pt 



y 



dy. 



It is not difficult to bound this integral away from 0, so it will remain to show that the integral is also bounded near 
0. Observe that for any yo > 0, 



yo P(ti < y) 



yo- sup 



dy - 



1 



0<j<yo y 



yo y 
P{zi <y) , 1 

Py'o 



dy 



The objective is now to establish that the term in braces can be bounded by a finite number By Theorem 2.5.2 in 
II25I . lim ^'^ ^ = 1, where f^:^ (y) is the density of Ti and 



' y^o Z{y) 



z{y) 



where A 



y' 



andC= (1 -j3)j3T^ > 0. 



To apply this result, we note that the derivative of the cumulative distribution function is the density, and apply 
L'Hospital's rule to obtain: 



= lim 



• lim 



ziy) 



lim ^^'^1 ^ lim 

y™+ y^+P y^+{l+li)yP y"0+ y"0+ (1+J3)/ 



7: lim 

1 + j3 y^0+ 



j2{l-P)' 



Therefore, there exists < 5 <yo such that for all < y < 5, < 1 and 



sup 

0<y<yo y 



P(Ti < y) 



i+P 



max sup 



P(Ti <y) 



v0<v<5 y 



l+p 



< max 1 , max 



1 



LetM = j3 



3'0(1 + 5tW)+7V 



Then for all «, E 



S<y<yoy^^ 



Pj^i <y) 

5<y<yo y^+P 
1 



< 1 



< fPM, as desired. 



□ 



Proof of Theorem^ We first show that (X"(f))„>i is a sequence of martingales. Note that the jumps have mean 
and finite first moment and for each n, E[Nnt] < °° for all f > 0. By the same techniques used in Lemma [T] we see 
that X"{t) is an ^"-martingale (and therefore a local martingale) for each n, where J^p is the natural filtration of X". 
Therefore, the first hypothesis of Proposition[T]holds. 
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Now we verify the second condition: sup„ E" 
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supj<, |AX"(,s')|] < oo for each t < °o. Fix f < °o. Observe that X"{s) 
is a cadlag step process with jumps of size -^^^jj at times G N. 

As in Theorem|2] the goal is to obtain for each f a uniform bound in n on the expectation. To this end, we first expand 
the expectation using the definition and standard technique of computing the integral of the tail probabilities. To apply 
Levy's Inequality, we symmetrize by replacing with ek^k where (e^.) is a sequence of i.i.d. Rademacher random 
variables, independent of {^i^). Then |<^<.| = \£k^k\ and e^^ii are i.i.d. in DOA{Z). Let S{m) denote the symmetrized 
sum. Since X"{s) has only finitely many jumps by time f, we can replace the sup up to time t by the max up to N„t and 
the following holds: 



£"[sup|AX"(^) 



m=0 



< 



P" { max >x] dx 

\<k<m cnP/'^ 
£k^k 



m=0 



2P" 
S{m 



k=i ' 



.p/2 



P{N„, - ; 



> x \ dx 



P{Nnt = m) 
P{Nnt = m) 



We can easily bound the sum for a finite number of terms. The goal is then to show that the tail of the sum is 

S{m) 



bounded. Since ekt,k G DOA{Z), then by 1 1 Exercise 9, p. 91], E 
S{m 



m > OTo, E 



(13) 



jl/2 



< /slZj + 1. We break up the sum at mo: 

2 mo oo 

Y^E"\S{m)\P{Nnt^m)+ £ E" 

m=0 m=mo+l 



,13/2 



J/2 



S{m) 



■ E\Z\, so there exists mq such that for all 



J/2 



m^l^P{Nnt = m) 



The sum up to mo is finite since E"-\S{m)\ < E\^i | < oo by the triangle inequahty. By applying the bound given 
above on the first moment of S{m)/m^l^, we can bound the tail of the sum by: 



(14) 



-(£|Z| + 1)£ 



1/2' 



Finally, we show Equation (fT4l) is finite. Since ^/x < x + 1 for all x > 0, and E [N„t /nP] < f^M by Lemma|2] we have 



(15) 



'nt 



We have shown that Equation ( fTST i is finite and is an upper bound for ^"[supj^, |AX"(5)|]. Therefore, the second 
hypothesis of Proposition[T]is verified. 

Since {X") satisfies UT, an application of Theorem 34 in lfT2l yields the conclusion that if {H",X") {H,B{E)) 
in (D](j2 , J I ) then there exists an appropriate filtration such that as n — ^ qq. 



H"{s-)dX'\s) 



f>0 



H{s-)dB{E,)\ in(DM,7i). 



□ 



The following example proves the conjecture of Germano et al. fS^: that the stochastic integral of the scaled 
compound Poisson process with respect to the scaled compound Poisson process converges weakly to the stochastic 
integral of Brownian motion with respect to Brownian motion. 

Example 2. Let (^,),gn be i.i.d. normal random variables with mean and variance 1, and (t,),^^ be i.i.d. exponential 
random variables with A = L ThenX(f ) = ^, is the usual compound Poisson process, and E{Nt) = t. Let X"(f ) = 



-'■=1 



be the scaled version. ThenX"(f) is a martingale for each n. Fix f < oo. We show ^"[sup^.^, |AX"(5)|] < o° 



using the same techniques as in the proof of Theorem [3] The sum in Equation dTSI) (with j3 = 1) up to mo can be 
bounded by a finite number and for m > mo, the tail is bounded by 



{E\Z\ + l)E 



1/2- 



< (£|Z| + l)Vf <oc. 



Since {X",X") ^ {B,B) in (Dr,/]) and the above shows the sequence {X") is good, then as n — )■ oo, 
(16) \ ['x"{s^)dX"{s)} =^ { f B{s) dB{s)] in(DR,7i). 

I Jo I ,>o I -'0 J 



r>0 



Remark 1. Meerschaert and Sheffler prove a general functional central limit theorem for CTRWs, in which the CTRW 
limit process is a time-changed Levy process A (£;) (see Theorem 4.2 in |15 |). This theorem holds in the Mi -topology 
but not in the stronger Ji -topology, except for the special case where A is Brownian motion (as discussed in Section 
|2]l. In proving Theorems |2] and |3] we rely on Kurtz and Protter's weak convergence theorem for stochastic integrals 
which holds in the 7i -topology. While we hope that a version of this theorem exists in the Mi -topology, this question 
remains open. 



4. Applications 

Knowing that a sequence of scaled CTRWs is good is useful in many situations. By definition, goodness of a 
sequence of scaled CTRWs {X") implies that if {H",X") converges jointly to {H,X), then the integral of H" driven 
by X" converges weakly to the integral of H driven by X. Since the stochastic integral driven by a scaled CTRW 
can be defined as a sum, the stochastic integral driven by a scaled CTRW can be easily simulated and used as an 
approximation for the stochastic integral driven by its limit process. Additionally, goodness is a necessary condition 
on the driving process for the weak convergence of stochastic differential equations to hold. Section|4T|focuses on two 
examples which apply the result of Theorem|3] in which the CTRW limit process is a time-changed Brownian motion, 
to cases of weak convergence of stochastic differential equations driven by scaled CTRWs. Section |42] describes the 
connection with particle tracking. 

4. 1 . Weak convergence of stochastic differential equations. The results in Section[3]can be combined with a theo- 
rem by Kurtz and Protter on weak convergence of solutions to stochastic differential equations lfT2l Theorem 38]. A 
key condition in this theorem is that the driving process is good. The conditions of the theorem are simplified in the 
case that the driving process is a scaled CTRW known to be good. We state an example in the case that the driving 
process satisfies the conditions in Theorem|3] where the hmit process is a time-changed Brownian motion. 

Example 3. Let {X")„>i be a sequence of semimartingales and (Z")„>i be a sequence of CTRWs satisfying the 
conditions in Theorem[3] Assume the following: 

• {X",Z") satisfies 

(17) X,"^ c/Z,". 

Jo 

• {X" , Z" ) is relatively compact in (Dr2 ,Ji). 

• are 7i -continuous. 

Then since the sequence (Z")„>i is good by Theorem[3]and Z" B{E) by Theorem[T] it follows from Theorem 38 in 
lfT2l that any limit point of the sequence {X") satisfies 

(18) Xt^ j'^F{X),^dBE(sy 



The next example involves a special kind of stochastic integral equation for which the form of the solution is known. 
Consider X, — \+ Jq X^- dY^, where F is a semimartingale and Yq =0. By 1 19_, Chapter 11, Theorem 37], the solution 
X is a semimartingale called the stochastic exponential of Y and is given by: 

(19) z, = exp|y,-i[}',y],| n (i+Ay,)exp|-Ay, + i(Ay,)2|, 

9 



where the infinite product converges. Using Theorem[3]and consequences of the goodness of the scaled CTRWs, we 
are able to verify directly, without checking any conditions on the stochastic integral equations, that the stochastic 
exponential of these scaled CTRWs converges weakly to the stochastic exponential of the CTRW limit process B{Et). 



N(nt) t 

Example 4. Let Z"(f) = Y where (^,) are i.i.d. Rademacher random variables and (t,) are i.i.d. strictly 

)3-stable random variables, j3 G (0, 1). Then Z" B{E) by Theorem[T]and (Z") is good by Theorem[3] Consider 

(20) X" = 1 + fx" dZ" 

Jo ' 

We want to verify that X" X in where X is a solution of the limiting equation 

(21) X,^l + f'x.^dBE,. 



We first analyze the form of the solutions. Since B{E,) is a semimartingale, the stochastic differential equation (ISTT l 
has a unique solution which is given by X, = exp{B£, ~ ^^^^^ ^" is a martingale, hence a semimartingale, 

with Z" (0) = 0, the stochastic exponential of Z" is the (unique) solution to ( l20b and is given by 

X"{t) = exp |z"(r) - ^[Z",Z"],| U (1 + AZ"(.))exp |-AZ"(.) + \{AZ"{s)f^ 
= W"it)-A"{t). 

To handle W"{t) — exp |z"(f) — i[Z",Z"],|, we appeal to the following two results; the first is an immediate 
corollary to Theorem 36 in 11121 and the second is Theorem 10.17 in |8|: 

1. If (Z") is good and Z"^Z in (Dr,7i), then (Z", [Z",Z"]) ^ (Z, [Z,Z]) in (ID)k2,/i). 

2. Let Y be an (^,) -semimartingale and (T,) be a finite continuous (^,) -time-change. Then [Y oT,Y oT] = [Y, Y] o T . 

Since Brownian motion B is a semimartingale, £ is a continuous nondecreasing process with E, <°o almost surely, 
and the quadratic variation [B,B], = t, then [B(E),B{E)], ^ {[B,B]oE), = E,. 

Since Z" ^ B{E) in (]I])r,7i), [Z",Z"] ^ [B{E),B{E)\ = £ in (©R,yi). Now using the 7i -continuity of addition on 
(Z", 5[Z",Z"]) and the continuous-mapping theorem, we have: 

Z" - i[Z«,Z"] ^ B{E) - i£ in {O^J,). 
By the 7i -continuity of the exponential and the continuous-mapping theorem. 



exp |z" - i[Z",Z"]| exp \b{E) - i/sj in (Pm,^i)- 



So W" ^ X in (Dk,-/i)- We need to show W" -A" ^X in (©r,7i). If A" ^ 1, then A" 1 in probability and so 
{A"',W") (1,X) in D X D. Since 1 and X are both continuous, multiplication here is continuous by 7] -continuity of 
multiplication ||23i . By the continuous-mapping theorem, W" A" ^ X. It remains to show A" 1, or equivalently 
logA" ^ 0. 
Observe that 

N(nt) , 

logA"(f)= £ [log(l+AZ"(.))-AZ"(.) + -(AZ"(.))2]. 

i=\ 

Also note that |AZ"(f)| = l-fyrl for t ~ and otherwise, since Z"(f) is a cadlag step process with jumps at times 
iM^keN. For n > 2^1^, \AZ"{s) 



1 1 

< —iTT^ < —, SO 



1 IA7"I^ 
|log(l+AZ"(.))-AZ"(.) + -(AZ"(.))2| < 

10 



by the alternating series estimation theorem. Therefore, 

N(nt) 



log(l+AZ"M)-AZ«(.) + -(AZ«(.))2 < ^ ^ < £ 
1=1 1=1 " 1=1 



Z'^'' 1 N{nt) 1 



„3i3/2 „/5 „/5/2- 



1 A^fnf) A^f«f) 1 

Since lim = and — £■(?) in {BuJi) as n ^ oo by Corollary 3.4 in yjj, — V^-r^ ^ in (Da,7i), 

which implies logA"(f) ^ 0, thereby confirming thatX" X. 

4.2. Particle tracking. These results also have applications to particle tracking. Particle tracking is a method of 
solving partial differential equations in cases where an analytic solution (i.e., closed form solution) cannot be found. 
Particle tracking involves first finding the stochastic differential equation (SDE) corresponding to the forward Kol- 
mogorov equation and second simulating the solution to the SDE ll24l . For example, a special case of Theorem 4. 1 in 
||6| considers the SDE: 

(22) dX, = a{X,)dE, +b{X,)dBE, M=x, 

and shows that the density p^{t^x,y) of Xt satisfies in the weak sense the following forward Kolmogorov equation: 



P it,x,y) 



(23) Dlp''{t,x,y) -- 

p''{0,x,y)^5,iy), 
where is the Caputo fractional derivative defined as 

Because of this correspondence between the SDE ( |22] | and the time-fractional differential equation (|23] l. particle 
tracking can be used to solve (l23l l. The solution to the SDE can be simulated using CTRWs where the waiting times 
arise from j3 -stable random variables as in Theorem[3] 
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